Formation of Ultracold Heteronuclear Dimers in Electric Fields 
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The formation of ultracold molecules via stimulated emission followed by a radiative deexcitation 
cascade in the presence of a static electric field is investigated. By analyzing the corresponding cross 
sections, we demonstrate the possibility to populate the lowest rotational excitations via photoasso- 
ciation. The modification of the radiative cascade due to the electric field leads to narrow rotational 
state distributions in the vibrational ground state. External fields might therefore represent an 
additional valuable tool towards the ultimate goal of quantum state preparation of molecules. 

PACS numbers: 33.80.Ps, 33.55.Be, 33.70.Ca 



Over the past decade investigations of ultracold quan- 
tum gases have been revealing a wealth of intriguing 
phenomena. In particular, ultracold molecular systems 
represent a paradigm including molecular Bose-Einstein 
condensates [2, 0, External fields are equally im- 
portant for the preparation and control of ultracold 
systems. They are used for cooling and trapping as 
well as for quantum state preparation or the tuning of 
atomic/molecular interactions, e.g., via Feshbach reso- 
nances [3, HI]. Moreover, they provide a tool to manip- 
ulate chemical reactions and collisions @, 0] , collisional 
spin relaxation in cold molecules and rovibrational 
spectra and transitions 0, [l(| • The quest for ultracold 
molecules being prepared in well-defined rovibrational 
quantum states is motivated by major perspectives such 
as the creation of quantum gases with novel many-body 
properties, ultracold state-to-state chemical reaction dy- 
namics, or molecular quantum computation (ill. Hil Il3| . 

A special focus is the study of heteronuclear polar 
dimers and the possibility to create dipolar quantum 
gases 14j. Due to the permanent electric dipole mo- 
ments of the constituent molecules, the intermolecu- 
lar interaction becomes long-ranged thereby introduc- 
ing exceptional properties for, e.g., the corresponding 
quantum gases. Electric fields play here an important 
role due to their immediate impact on polar systems 
which leads to rotational orientation and alignment of 
the molecules. The control of the formation process of 
ultracold polar molecules is a key issue to arrive at molec- 
ular quantum gases. Photoassociation of ultracold atoms 
is a widespread technique to produce ultracold molecules 
15 1 . Indeed, the formation via photoassociation of KRb 
3, LiCs NaCs [3 and RbCs [H, [U has been 
reported recently. 

Let us consider a mixture of two atomic species in their 
electronic ground states, exposed to a linearly polarized 
laser with a frequency corresponding to the transition 
from the continuum to a highly excited vibrational bound 
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FIG. 1: (Color online) Electronic potential energy curve and 
dipole moment function of the electronic ground state X E + 
of the LiCs molecule, in atomic units. The schematics of the 
OPA and RDC are also shown (not to scale). 



state of the electronic ground state. The corresponding 
process of stimulated emission of a single photon allows 
for the formation of a polar molecule, see fig. 1 (note that 
hyperfine states are not considered) . It has been recently 
studied for the LiH and NaH molecules 21 , [22| focusing 
on the ultracold regime where s — > p-wave transitions 



dominate and only rovibrational bound states with an- 
gular momentum equal to one can be produced. Subse- 
quently, a radiative cascade of rovibrational transitions 
leads to the vibrational ground state, thereby populating 
a broad range of rotational molecular states. 

In the present work we explore the impact of an addi- 
tional homogeneous static electric field on the formation 
process of polar dimers. As indicated above, the latter 
consists of the one-photon stimulated association (OPA) 
followed by a radiative deexcitation cascade (RDC). We 
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demonstrate that the presence of the static field allows 
one to populate via the OPA electrically dressed states 
evolving from field-free levels with zero angular momen- 
tum. This in combination with the following RDC yields 
a final rotational state distribution within the lowest vi- 
brational band which is significantly narrower compared 
to the field-free case. Consequently, a static electric field 
can serve as a tool to prepare a molecular gas with only a 
limited number of well-defined molecular quantum states 
being populated. This might be a helpful step towards 
the goal of reaching quantum gases of molecules in their 
rovibrational ground states. 

Let us assume that perturbation theory suffices for the 
description of the interaction with the laser field and for 
the interaction of the electronic motion with the static 
electric field. However, a nonperturbative treatment is 
indispensable for the impact of the electric field on the 
nuclear dynamics. The Hamiltonian for the rovibrational 
motion of the diatomic system in the Born-Oppcnhcimcr 
approximation reads 
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-e(R)—FD(R) cos 9 
(1) 

where the molecule fixed frame with the origin at the 
center of mass of the nuclei has been employed. (R, 6, <j>) 
are the internuclear distance and the Euler angles, re- 
spectively, /it, J(9,4>), e(R), F, and D(R) are the re- 
duced mass of the nuclei, the rotational angular momen- 
tum, the field-free electronic potential energy curve, the 
electric field strength, and the electronic dipole moment 
function, respectively. The electric field is parallel to 
the z-axis of the laboratory frame. In field-free space, 
each bound state of the molecule is characterized by its 
vibrational, rotational, and magnetic quantum numbers 
(v,J,M). In the presence of the electric field only the 
magnetic quantum number M is conserved. However, for 
reasons of addressability we will refer to the electrically 
dressed states by means of the corresponding field-free 
quantum numbers. 

In the framework of the dipole approximation, the 
cross section of the stimulated emission process from the 
continuum \&£(R) to the rovibrational state ^(R) is 
given by 



tt(E - E a ) 
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where a labels the bound state and E, E a are the contin- 
uum and rovibrational energies, respectively. The con- 
tinuum wavefunctions are energy normalized while the 
rovibrational ones are C 2 normalized. 

The nuclear equation of motion associated with the 
Hamiltonian |T]) is solved by means of a hybrid com- 
putational approach together with a Krylov-type diag- 
onalization technique [23j. For the vibrational coordi- 
nate we use a mapped discrete variable representation 



based on sine functions 24] in a box of size L. A ba- 
sis set expansion in terms of the associated Legendre 
functions is performed for the angular coordinate. This 
computational technique provides an accurate descrip- 
tion of highly excited bound states, while the continuum 
spectrum is discretized and described by C 2 normalized 
states. We emphasize that the traditional technique to 
represent energy normalized continuum wave functions 
by means of C 2 normalized functions 25j is here not ap- 
plicable due to the field-induced coupling between the ro- 
tational and vibrational motions. Moreover, the coupling 
significantly alters the definition of the density of contin- 
uum states. This issue has been solved by employing 
the time-dependent formalism. Since our study focuses 
on the ultracold regime with temperatures below 1 mK, 
very large discretization boxes are required to properly 
describe the continuum. 

Being a topical example, we focus on the polar alkali 
dimer LiCs. The potential energy curve of its X 1 E + 
electronic ground state is taken from experimental data 
[26| and the electric dipole moment function from semi- 
empirical calculations 27J (see fig. 1). Continuum ener- 
gies corresponding to T = 10, 100, 500 fiK are considered. 
We analyze the stimulated emission process for the pop- 
ulation of the rovibrational levels (44, J, 0), which is a 
reasonable but robust choice, with J = 0,1, requiring 
laser wavelengths of 181.13 /zm and 181.56 /im, respec- 
tively. Similar results are obtained for neighbouring vi- 
brational bands. The dependence of the total cross sec- 
tions on the electric field strength is shown in fig. 2. It is 
worth noting that the considered regime of field strengths 
covers the experimentally accessible range and beyond: 
Static fields F > 4 • 10~ 5 a.u. = 200 kV/cm are difficult 
to achieve in the laboratory. The cross sections increase 
with increasing temperature due to a larger overlap be- 
tween the continuum and bound states. In the absence 
of the electric field, the ultracold regime is dominated by 
s —> p-w&ve transitions and the cross section for popu- 
lating the (44, 0, 0) state via the OPA process is several 
orders of magnitude smaller than the corresponding cross 
section involving the (44, 1,0) state. 

Augmenting the field strength, the cross sections 
change significantly which can be explained by the hy- 
bridization of the angular motion. For F > 3 ■ 10~ 6 a.u., 
the cross section for the (44, 0, 0) state becomes larger 
than the one belonging to the (44,1,0) state. Adopt- 
ing T = 10 fxK, the cross section of the J = state in- 
creases monotonically by two orders of magnitude within 
the regime F = 10~ 7 — 10~ 5 a.u. Further increasing the 
field strength, it exhibits a plateau followed by a weakly 
pronounced minimum and a strong increase thereafter. 
In contrast to this, the cross section of the (44, 1, 0) state 
shows a plateau in the weak field regime. For stronger 
fields it rapidly decreases and reaches a deep minimum 
for F = 1.1 • 10~ 5 a.u., followed by a broad maximum and 
a significant decrease for strong fields. This behaviour of 
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FIG. 2: (Color online) Stimulated emission cross sections as a 
function of the electric field strength for the states emerging 
from field-free levels with quantum numbers (44, 0, 0) (dotted, 
dotted-dashed and double dotted-dashed lines) and (44, 1,0) 
(solid, dashed and short dashed lines) for continuum energies 
corresponding to T — 500, 100 and 10 [iK, respectively. 



a laser beam with intensity 1 ^X-, and a temperature 
T = 1 mK we obtain a formation rate of 10 4 molecules 
per second. Recently, the formation of LiCs molecules in 
a two-species magneto-optical trap by a two-photon pro- 
cess has been reported [17j. The molecules are formed 
in the electronic ground state, but an analysis of the fi- 
nal vibrational state distribution is not provided. With 
densities of 10 10 cm" 3 and 5 • 10 9 cm" 3 for the Li and 
Cs atoms, respectively, and assuming a temperature of 
100 fiK, a molecular production rate between 1.4 ± 0.8 
and 140 ± 80 molecules per second is estimated. Using 
these densities for the Li and Cs atoms, our theoretical 
results for the molecular formation rate via the OPA pro- 
cess are in the same order of magnitude. The reverse pro- 
cess of one photon absorption leading to the dissociation 
of the molecules is suppressed by the RDC (see below) , 
although this is much more pronounced for light hydrides 



21] compared to alkali dimers. It might be further elim- 



inated by applying (chirped) laser pulses. Effects due to 
vibrational quenching can be neglected for not too large 
molecular densities. 

Since the photoassociation process results in vibra- 
tionally highly excited states, a cascade of spontaneous 
emission processes will follow. The overall transition 
probability per unit time reads T a = Y^vU^f'" ^a,a' 
where the summation includes all open decay channels, 
i.e., final levels a' with M' € {M,M ± 1}. The corre- 
sponding transition rates are 

Ta ia > = ^^\{* a \D{R)f M ,{6)\* a ,)\* (3) 

with fAi(0) — cos 9 and /m±i(^) — sin#. hui a a r is 
the energy difference between the initial and final states. 
The total and single-channel radiative lifetimes are r Q = 
(L Q ) _1 and T a<a > — (r Q , Q ') _1 ; respectively. Our aim is 
to analyze the final rotational state distribution in the 
lowest vibrational band resulting from the RDC. Start- 
ing from an initial state we therefore add up all decay 
paths resulting in a specific final state. This provides us 
the final J, M-state distributions. 

Besides the possibility to form molecules evolving from 
field- free states with zero angular momentum by applying 
an electric field, the field interaction significantly alters 
the radiative properties of the bound molecular states. 
The rotational state distribution with varying (J, M) in 
the ground vibrational state which results from the ra- 
diation cascade of the (44,1,0) state in the absence of 
the field is illustrated in fig. 3. The inset shows the cu- 
mulative population for fixed rotational quantum num- 
bers J (\M\ < J). Analogous results for the initial state 
ao = (44, 0, 0), but for F = 4- 10~ 5 a.u., are presented in 
fig. 4. It should be noted that negative magnetic quan- 
tum numbers obey the same distribution due to the de- 
generacy of ±M states and the initial value M — 0. Aug- 
menting the field from zero to F — 4 ■ 10~ 5 a.u., corre- 
10 12 cm -3 , a volume V — 10~ 6 cm 3 illuminated by sponding to 200—, the radiative lifetime of the (44,0,0) 



the cross sections is due to the dominance of s — > p- 
wave transitions in the ultracold regime combined with 
the hybridization of the angular motion of the rovibra- 
tional state. Indeed, the presence of the plateaus in both 
cross sections is accompanied by an analogous constancy 
of the contribution of the p-wave to the corresponding 
wave functions. Similarly, the minimum for the J = 1 
state at F = 1.1 • 10" 5 a.u. is due to the dominance of s, 
d, and /-waves along with a rather weak contribution of 
the p-wave. 

Let us shortly analyze the hybridization of the angular 
motion for these two states in more detail. The (44,0,0) 
state has a high-field-seeking character with an increas- 
ing expectation value (J 2 ) for augmenting fields, i.e., the 
contribution of higher field- free rotational states becomes 
larger as the electric field strength is enhanced. Conse- 
quently, this level shows a pronounced angular momen- 
tum hybridization J 2 (44, 0, 0) = (J 2 )(4 4 ,o.o) - J(J+l) = 
2.08 and 7.88, for F = 10" 5 a.u. and KT 4 a.u., respec- 
tively, with J being the corresponding field-free rota- 
tional quantum number. On the other hand, the (44, 1, 0) 
state shows initially a low-field-seeking behaviour where 
the admixing of lower rotations dominates and (J 2 ) de- 
creases: J 2 , (44, 1,0) = -0.11 for F = 10" 6 a.u. Further 
enhancing the electric field, this level becomes a high- 
field-seeker, and (J 2 ) increases after reaching a minimum. 
In the strong field regime, this level exhibits a very pro- 
nounced hybridization J 2 (44, 1, 0) = 21.45 for F = 10" 4 
a.u. 

The formation rates in the OPA process can be varied 
significantly by simply changing the intensity of the ap- 
plied laser. Specifically, if we assume an atomic density 
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FIG. 3: (Color online) Population distribution as a function 
of the magnetic quantum number M and the difference J — M 
for the decay cascade starting from the (44, 1, 0) state in the 
absence of the electric field. The inset shows the cumulative 
population for fixed rotational quantum numbers. 



state increases from 1.57 s to 2.15 s, while for the (44, 1, 0) 
state it is reduced from 1.87 to 1.76 s. In both cases, the 
prevailing variation of the vibrational quantum number 
for each transition changes throughout the cascade from 
Av = 6 for the v = 44 state to Av = 1 for low-lying vi- 
brational states. It is important to note that, in contrast 
to its favourable large electric dipole moment, LiCs rep- 
resents a molecule with quite long radiative decay times. 
In comparison, heavier alkali dimers have even longer life- 
times, e.g., kiloseconds for the KRb molecule [28j while 
light hydrides such as LiH exhibit much shorter lifetimes 
of the order of a few milliseconds [2l[ . 

In the absence of the field, the selection rules A J = ±1, 
AM £ {0, ±1} for dipole transitions hold, giving rise to a 
very broad distribution of the population in the lowest vi- 
brational band: a large number of states exhibit a similar 
population. For a fixed J, the fully angular-momentum- 
polarized states always possess the largest population, 
and for fixed M the population slowly decreases with in- 
creasing J. Moreover, the two fastest paths populate the 
rotationally highly excited (0, 20, ±20) states in 159 s. 
The cumulative population of the rotational bands ex- 
hibits a broad distribution with a maximum at J = 5. 
Similar results are obtained for the (44, 0, 0) state. Once 
the vibrational ground state is populated, the rotational 
cascade is an extremely slow process: The lifetimes of the 
(0, 1,±1) and (0,5, ±5) levels are 1.88 • 10 6 s and 1.1 ■ 10 4 
s, respectively. We remark that heating due to the ran- 
dom direction of the emitted photons in the course of 
the RDC is a minor effect amounting to several tens of 
nanokelvins for a typical setup [2l| . 

The population of the final rotational states in the 
lowest vibrational band following the OPA and RDC 
changes drastically in the presence of the electric field. 
Due to the hybridization of the angular motion, only 
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FIG. 4: (Color online) Same as in fig. [3] but with (44, 0, 0) as 
initial state and for F = 4 ■ 10~ 5 a.u. 



the selection rule AM £ {0, ±1} for the magnetic quan- 
tum number holds and many new transitions are possi- 
ble. In particular, for fully polarized angular momen- 
tum states the importance of purely vibrational tran- 
sitions AJ = AM = increases with increasing field 
strength and these transitions become dominant in the 
strong field regime [Io| . Consequently, the distribution 
is much narrower and exhibits pronounced peaks for a 
few fully angular-momcntum-polarized states as can be 
seen by comparing figs. 3 and 4. For F = 4 • 10~ 5 
a.u., the RDC is dominated by these purely vibrational 
transitions and the rovibrational ground state possesses 
the largest population of 5.75% compared to 2.02% for 
F = 0. This effect can be further enhanced by employing 
even higher fields. In addition, the rovibrational ground 
state is populated by the most probable paths. Accord- 
ingly, we observe in the inset of fig. 4 that the maximum 
of the cumulative rotational state population is shifted to 
the J = 2 band and includes 17.4% of the overall norm. 
The entire RDC process to the vibrational ground state 
typically takes several minutes. We remark that the sub- 
sequent rotational decay itself is enhanced significantly 
in the presence of the field and the corresponding fully 
angular-momentum-polarized states possess a unique ra- 
diative decay route to the rovibrational ground state [lfj| . 

The above analysis demonstrates that the formation 
of ultracold molecules in their electronic ground state 
via a single-photon stimulated emission process followed 
by a radiative cascade can be significantly altered by the 
presence of a static electric field. The resulting narrow 
rotational distributions in the lowest vibrational band 
represent a remarkable step towards the ultimate goal 
of degenerate quantum gases of molecules in their rovi- 
brational ground state, or more generally in any desired 
unique quantum state. Although our study focuses on 
the LiCs molecule, a qualitatively similar behaviour has 
to be expected for other heteronuclear dimers depending 
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on the corresponding potential energy curve and electric 
dipole moment function. 
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